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Abstract 



We measure the CP-violating asymmetries in decays to the D" -^ K^K^ and D" ^ tt+tt^ CP eigenstates using 
540 fb~ of data collected with the Belle detector at or near the T(4S') resonance. Cabibbo-favored D —^ K~tv'^ 
decays are used to correct for systematic detector effects. The results, Aqp = (—0.43 ± 0.30 ± 0.11)% and AJp = 
(+0.43 ± 0.52 ± 0.12)%, are consistent with no CP violation. 
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1. Introduction 

Decays of neutral D mesons are a promising area 
in which to search for physics beyond the Standard 
Model (SM). Recently, evidence for mixing in this 
system has been obtained |1I2I3J . However, whether 
the effect observed is due to the Cabibbo-Kobayashi- 
Maskawa (CKM) theory or due to new physics (NP) 
has yet to be determined and will require further 



measurements to resolve. One possible measurement 
sensitive to NP is that of a CP asymmetry in D° 
decays to Cabibbo-suppressed (CS) final states [4]. 
Within the SM such an asymmetry is predicted to 
be very small (< 0.1%), but within NP scenarios it 
can be substantial (> 1%) [415] . 

In this Letter we present a high statistics search 
for a CP asymmetry in the CS modes D° — > K~^K~ 
and £)° -^ tt+tt". These final states are accessible 



to both D'^ and Z)" mesons. The time-integrated 
CP asymmetry for decays into a CP eigenstate / is 
defined as 






/) - r(5" ^ /) 



/) + V{D^ ^ /) 






+ aind • (1) 



This quantity receives contributions from both di- 
rect (a;5) and indirect (oind) CP violation (CPV) 4 . 
While the direct contribution is in general distinct 
for different final states, the indirect contribu- 
tion is the same. The indirect CPV contribution 
is constrained by our recent measurement of the 
lifetime difference using D^{D^) — * K^K^ ,Tr^Tr^ 
decays [1]: Ar = -a-.nd = (0.01 ± 0.30 ± 0.15)%. 
CP asymmetries in CS decays have been searched 
for previously using D° -^ K^K^ ,n^TT~ [5] and 



D" 



Tr+Tr-Tr°,K+K- 
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2. Method 



The flavour of neutral D mesons at production 
is tagged by reconstructing Z?*+ -^ D'^'^t "^c- 
cayfUJ in which the charge of the low momentum 
pion, TTs, determines the flavour of the D^ me- 
son. The measured asymmetry, A{^^ — [N{D'^ -^ 
f)-N{D^-. f)]/[N{D^ ^ f) + NiD° ^ /)], with 
/ = K'^K^ ,Tr^Tr~ and N denoting the number of 
reconstructed decays, can be written as a sum of 
several (assumed small) contributions: 



aL = Afr + aI.^ + a: 



CP 



(2) 



In addition to the intrinsic asymmetry A^p there 
is a contribution due to an asymmetry in the re- 
construction efRciencies of oppositely charged tt^ 
{A'^). Since the final state / is self-conjugate, its 
reconstruction efficiency does not affect A(^^. Fur- 
thermore, there is a forward-backward asymmetry 
{Aps) in the production oi D*^ mesons in e^e^ —>■ 
cc arising from 7 — Z" interference and higher order 
QED effects [5]. This term is an odd function of 
the cosine of the £>*+ production polar angle in the 
center-of-mass (CM) systeno (cos^*). Since our 
detector acceptance is not symmetric with respect 
to cos 9* , the measurement is performed in bins of 
cos 9* . This allows us to correct for acceptance and 
extract both ApB and A^p as described below. 



^ Charge conjugated processes are implied throughout the 
paper, unless explicitly noted otherwise. 
^ Symbols with an asterix in the paper denote quantities in 
the CM frame, while those without asterix denote quantities 
in the laboratory frame. 



To reliably determine A'^ we adopt the method of 
Ref. 6 with some appropriate modifications. In ad- 
dition to the D° -^ h^h^ modes mentioned above, 
we also reconstruct two 1)° -^ K^tt^ samples: one 
consisting of D mesons with tagged initial flavour, 
and one consisting of untagged candidates. The mea- 
sured asymmetries for these modes can be written 
as 
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CP 
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FB + ^CP + ^£ ■ (3) 

A notable difference with ([2]) is that this final state is 
not self-conjugate and thus an additional term A^'^ 
appears as a consequence of a possible asymmetry 
in the reconstruction efficiency. We first use the two 
measurements in ^ to determine A^; we then in- 
sert the result into ^ and use the fact that Aps 
is antisymmetric with respect to cos 9* and A^jp is 
independent of this variable. 

Reconstruction efficiencies and their asymme- 
tries A*, however, are functions of momenta of 
particles z = tt^. Kit in the laboratory frame. For 
a D^ meson with a given momentum pjja , the ef- 
ficiency of reconstructing the final state isT" 7r+ is 
ex-KiPDo) = J eK{pK)e^{PiT)wj;^g{pK,P7r)dpKdp-„, 
where w^^g {pk , P-k) denotes the 6-dimensional dis- 
tribution of final state particles. For a given puo, 
this distribution is independent of whether the D 
meson candidate is flavour-tagged or not. Using the 
same selection criteria for the K and tt candidates 
in the tagged and untagged sample imposes equality 
of the selection efRciencies £K(TT){PK(Tr))- Hence the 
asymmetry A^'^{pua) is identical for tagged and 
untagged D mesons of a given momentum pijo , as 
implied by ([3|) . Since the distribution of D° mesons 
is uniform in the azimuthal angle the dimension of 
the problem can be reduced. It is sufficient to obtain 
A^'^ as a function of the magnitude and polar angle 
of the laboratory momentum, pijo and cos^do. 

The slow pion asymmetry A^ depends on its mo- 
mentum p^^ and is independent of the D° final state. 
Since the tTs azimuthal angle distribution is also 
found to be uniform, A'^ is examined as a function 
of (p^^, cos 6*^ J. 

3. Measurement 

The measurement is based on 540 fb^"'^ of data 
recorded by the Belle detector [5] at the KEKB 
asymmetric-energy e~^e~ collider [lOj . running at 
the CM energy of the T(4iS') resonance and 60 MeV 



below. The Belle detector is described in detail else- 
where |9J: it includes in particular a silicon vertex 
detector (SVD), a central drift chamber, an array 
of aerogel Cherenkov counters, and time-of-flight 
scintillation counters. Two different SVD configura- 
tions were used: a 3-layer configuration for the first 
f53 fb~^ of data, and a 4-layer configuration [TT] for 
the remaining data. 

We reconstruct D*+ -^ D°tt+, D° ^ K+ K- , 
K^t:^, tt+tt^ decay chains, as well as the decay 
D" -^ K~TT^ without requiring an accompanying 
D*^ decay. Each final state charged particle is re- 
quired to have at least two associated SVD hits in 
each of the two measuring coordinates. To select 
pion and kaon candidates, we impose standard par- 
ticle identification criteria [12]. Z?" daughter parti- 
cles are refitted to a common vertex. The D^ pro- 
duction vertex is found by constraining the D'^ (and 
TTs for the tagged decays) to originate from the e+e^ 
interaction region. Confidence levels exceeding 10^"^ 
are required for both fits. The Z)*+ (Z?*^ for untagged 
decays) momentum must satisfy pjj > 2.5 GeV/c^ 
in order to reject Z3-mesons produced in _B-meson 
decays and to suppress combinatorial background. 

We accept candidates with a Z?" invariant mass 
M in the range 1.81 GeV/c^ < M < 1.91 GeV/c^ 
For final states with a tt^ , we require that the energy 
released in the £)*+ decay, q — [Mjj,+ — M — m.,^)c^ , 
be less than 20 MeV. In this expression, Mjj*+ is 
the invariant mass of the D^n^ combination and m^ 
is the charged pion mass. For the small fraction of 
events with multiple candidates (0.1% for the tagged 
samples, 2.9% for the untagged sample), we select 
only one candidate: that in which the sum of the pro- 
duction and decay vertex x^'s is smallest. We also 
requirg ^ 1 1 cos0£)o| < 0-9 to remove events in which 
large slow pion asymmetry corrections and conse- 
quently large systematic uncertainties are expected. 
The resulting invariant mass spectra are shown in 
Fig. [11 

We measure the signal yield by performing a 
mass-sideband subtraction, as this method is ro- 
bust and reduces sensitivity to the signal shape. 
The possibility of a non-linear background shape is 
considered as a systematic uncertainty. The sizes 
of signal windows in M and q are chosen to min- 
imize the expected statistical error on the Acp 
measurement. Using the Monte Carlo (MC) sim- 
ulation, which has been tuned to reproduce the 
signal shapes and the signal-to-background ratios 
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Fig. 1. Invariant mass spectra of selected events. For the 
tagged data samples (b,c,d) events with [Aql <1 MeV are 
selected. The cross-hatched area represents the signal region; 
the sideband positions are indicated by vertical lines. 

of the real data, the optimal signal windows are 
found to be |AM| < 17.3 (18.6, 16.8) MeV/c^ and 
|Ag| < 1.00 (1.85,0.90) MeV for the KK (Xtt, tttt) 
final states. The quantities AM and Aq measure the 
difference of the corresponding observable and the 
nominal D'^ mass and the nominal released energy 
of the D*'^ decay, respectively. Sidebands of the 
same size as signal window are chosen starting at 
±20 MeV/c^ from the D'^ nominal mass. Within the 
optimal signal window we find 6.3 x 10^ untagged 
K^TT^ signal events with a purity of 80%; the num- 
ber of tagged signal events is 120 x 10'^ K~^K~, 
1.3 X 10*^ K-Tr+ and 51 x 10^ tt+tt", with purities 
of 97%, 99% and 91%, respectively. 

We determine first the asymmetry A™*^^ of the 
untagged Ktt sample in 20 x 20 bins of the two- 
dimensional phase space {pjjo , cos 0]jo ) by 



A. 



untag ^^^J 



N,. 



N,, 



N,.. 



(4) 



This cut limits the range of measurement to | cos 6* \ < 0.8 



where Nij and Nij are the numbers of reconstructed 
D^ and D^ decays, respectively, in bin ij. In order 
to avoid large statistical fluctuations near the phase 
space boundaries, we calculate the asymmetry only 
for those bins having Nij + Nij > 1000. This asym- 
metry is used to correct the tagged Kn events by 
weighting each D'^{D'~') candidate falling into a valid 
bin with a weight 
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Fig. 2. Asymmetry of the slow pion efficiency, A^, in mo- 
mentum slices for the 3-layer (left) and the 4-layer (right) 
SVD configurations. 
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Other candidates are discarded. The weighting ap- 
phed to the tagged Ktt decays results in a measured 
Ajpf free of aU contributions in ([3]) except for A^ . 

The slow pion asymmetry in bin kl of the phase 
space {ptt^ , cos 6*^^ ) is thus determined with 



A^ — 



nki - riki 
nu + flki 



(6) 



where n^; (jiki) are the sums of weights of the D^{D^) 
candidates falling in that bin. Again, we consider 
only bins with rifci -|-nfci > 1000. The resulting asym- 
metry A^ (ptt^ , cos Ott^ ) determined in 5 x 5 bins for 
the two SVD configurations is shown in Fig. [21 Aver- 
aging over the phase space the correction due to the 
slow pion efficiency is found to be (-1-0.76 ± 0.09)%. 
The slow pion asymmetry is used to correct the 
KK and tttt events. The D'^ / D^ candidates are 
weighted according to 



wd° = l~ A'^ (p^^ , cos 0^^ ) , 

W^O ^l+A'^ {Py,^ , cos Oy,^ ) , 



(7) 



and only candidates in bins with valid A'^ measure- 
ments are taken into account. This procedure results 



in a corrected asymmetry A(^^ of ^ , A(^l°" , which 
is free of the contribution due to the slow pion effi- 
ciency asymmetry. It is calculated as 



Ai,r-(cosr) = -:;--:; --;;----:; , (s) 



mJ^ (cos 9*) + m' (cos 0* ) 



where ra^ (rn^ ) represent the sum of weights of the 
D^{D'^) candidates in each bin of cos 6**. 

Finally, taking into account their specific depen- 
dence on cos 0* , the asymmetries Acp and Aps are 
extracted by adding or subtracting bins at ± cos 9* : 



A^ 



A(ff"{cos9*) + Aff^°"{~ cos 9*) 



ApB — 



A/'™"(cosr)-^/<;™"( 



COS! 



') 



(9) 



The results are presented in Fig. [3l By fitting a 
constant to the A^p data points we obtain results 
consistent with no CP violation: 



A^P =(-0.43 ±0.30)% 
(+0.43 ±0.52)% 



^CP 



(10) 



The errors are statistical only; however, the statisti- 
cal uncertainties of the slow pion corrections are not 
included. The forward-backward asymmetry ApB 
decreases with cos 6** and has a value « —3% at 
cos 61* =0.8; results from the two samples are con- 
sistent. At leading order, the asymmetry at this en- 
ergy is expected to be yl^^ (cos 0*) = a'^'^cos0*/(l + 
cos^ 9*), with a"^ = -2.9% [13]. A simultaneous fit 
to the two samples yields an acceptable goodness- 
of-fit (xV^dof = 4.5/7) and a'=^ = (-4.9 ± 0.8)%, 
where the error is statistical (see Fig.|3|). 

4. Systematics 

The experimental procedure was checked using 
the generic continuum MC simulation; the resulting 
Acp and App were found to be in a good agree- 
ment with the generated values. We also tested for 
possible bias in the result by re-weighting MC sam- 
ples with several non-zero Acp values; no significant 
bias was found. 

We consider three sources of systematic uncer- 
tainty to be significant (Table [1]) . The first source 
is the mass-sideband subtraction procedure used for 
signal counting. Possible systematic uncertainties 
arise due to the difference in signal shapes of I?" 
and Z)" candidates and due to the possible differ- 
ence in the background between the signal window 
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Table 1 

Summary of systematic uncertainties in Acp. 



Fig. 3. CP-violating asymmetries in (a) KK and (b) tttt fi- 
nal states, and forward-backward asymmetries in (c) KK 
and (d) tttt final states. The solid curves represent the cen- 
tral values obtained from the least square minimizations; the 
dashed curves in (c) and (d) show the leading order expec- 
tation. 

and sideband. The former source can introduce an 
additional asymmetry if the signal window is not 
sufRciently wide. We observe small but significant 
differences in the q signal shape of the tagged sam- 
ples. By studying the normalized (in order to asses 
only the effect of the shape difference) q distribu- 
tions of the tagged D^{D^) -^ Kir samples we esti- 
mate the systematic uncertainty of this source to be 
0.02% (0.04%) for the KK (tttt) sample. To account 
for a possible difference in backgrounds we vary the 
position of the sideband. We find 0.01% {KK) and 
0.03% (tttt) variations in the result. Background due 
to a correctly reconstructed D^ candidate combined 
with a random slow pion is not removed by the 
M sideband subtraction. Its fraction (0.6%) is esti- 
mated from the tuned MC simulation. The possible 
asymmetry induced by this type of background is 
estimated from the q sideband to be at most 0.03%. 
The second source of systematic error is the slow 
pion efficiency correction. The statistical errors on 
-^^{Pt^s 1 cos Ot^^ ) contribute an uncertainty of 0.09%. 
The impact of binning of the slow pion asymmetry 
is studied by producing maps with three different 
choices of bin sizes (10 x 10, 20 x 20, 50 x 50 for 
A™*^s, and 5x5, 10 x 10, 20 x 20 for A^) and 
repeating the procedure for extracting Acp ■ We find 
0.03% (KK) and 0.02% (tttt) variations in the result. 



Source 



D° -^ K+K- D" 



■ TT ' TT 



Signal counting 
Slow pion corrections 
Aqp extraction 



0.04% 
0.10% 
0.03% 



0.06% 
0.10% 
0.04% 



Sum in quadrature 



0.11% 



0.12% 



The minimum required number of events per bin is 
varied from 100 to 10000, and the resulting variation 
in Acp is 0.04% (0.03%) for the KK (tttt) sample. 

The third source of systematic uncertainty is the 
Aqp extraction procedure. By varying the binning 
in I cos 0*1 we obtain a 0.03% variation in the re- 
sult. We change the treatment of the running peri- 
ods with 3- and 4-layer SVD configuration; we find 
an 0.01% (0.02%) change in the resuh for the KK 
(tttt) sample. 

Finally, we add the individual contributions in 
quadrature to obtain the total systematic uncer- 
tainty. The result is 0.11% (0.12%) for the KK (tttt) 
sample. The dominant source is the statistical un- 
certainty on A'^ , and thus the majority of the sys- 
tematic error will decrease when a larger Kt: data 
sample is available. 

5. Conclusions 

We measure time-integrated CP-violating asym- 
metries Acp in decays to CP eigenstates D^ -^ 
K+K- and D° -^ tt+tt" using 540 fb'^ of data. 
The detector-induced asymmetries are corrected 
with a precision of 0.1% by using tagged and un- 
tagged D° -^ K~n+ decays. We obtain: 

tKK 



A^P = (-0.43 ± 0.30 ± 0.11)% , 



A TTTT ^ 

aKK 
^CP 



(+0.43 ±0.52 ±0.12)% , 

-0.86 ±0.60 ±0.07)%. 



AITTT 

^CP 



(11) 



The results show no evidence for CP violation and 
agree with SM predictions. In (fTTjl we also list the 
difference A^p — A'^J'p, which is calculated by treat- 
ing the systematic errors arising from the slow pion 
corrections and Acp extraction as fully correlated 
between the two modes. A significant difference be- 
tween the measured asymmetries in the KK and tttt 
modes would be a sign of direct CPV (Eq. ([T])). 

To determine the direct CPV asymmetries aj^ of 
^ , the results in (|lip can be compared to the result 
for the indirect CPV asymmetry in Rcf. Ij. While 
the selected data samples of D'-* -^ K+K^ ,ti+tt~ 



in the two measurenients are almost identical, the 
methods of extracting the CP violating asymme- 
tries depend on different observables and hence the 
statistical uncertainties are uncorrelated. The same 
holds also for the systematic errors. The direct CPV 
asymmetries following from the sum of A^p and Ar 
are 



JiK 



(-0.42 ±0.42 ±0.19)%, 
= (+0.44 ±0.60 ±0.19)% . 



The measurement uncertainties are above the level 
of the expected asymmetry in the SM. 

We also measure the forward-backward asymme- 
try in the production of D*+ that arises from the 
underlying asymmetry in the e^e'' —>■ cc process. 
The asymmetry agrees with the form A'^g (cos 6*) — 
a'^'^ cos 6* / {1 + cos^ 6*) expected at leading order, 
but we find a"^^ = (-4.9 ± 0.8)%, larger than the 
leading-order value of —2.9%. Radiative and other 
(hadronic) corrections are expected to cause the ef- 
fective a^^ to deviate from its leading-order value. 
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